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Abstract : Several theoretical models of earth's albedo radiation have been pro- 
posed by different authors. By comparing disturbing accelerations computed from 
a model to those measured in flight with the CACTUS Accelerometer, we may exa- 
mine its likelihood and modify it accordingly the results. 

Computation of the satellite orbit perturbations from a model is very long because 
for each position of this satellite we have to sura up the fluxes coming from each 
elementary surface of the terrestrial portion visible from the satellite. We in- 
crease the speed of computation ten times without significant lost of accuracy 
thanks to a stocking of some intermediate results. Now it is possible to confront 
the orbit perturbations computed from the selected model with the measurements 
of these perturbations found with satellite as LAGEOS. 

Conclusions are produced. 
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1. Choix d'un modele de rayonnement d ’albedo terrestre 


L'etude des accelerations mesurees par 1 ' accelerometre CACTUS embarque a bord du 
satellite CASTOR a permis de mettre en evidence et d'analyser la pression de ra- 
diation (et tout particulierement sa composante verticale) ayant pour origine le 
rayonnement infra-rouge emis par la Terre, et les rayonnements solaires rediffuse 
et reflechi par celle-ci. (BERNARD A. et al. 1977, LALA P. 1978, WALCH J.J. 1981). 

II devenait alors possible d'adopter un modele pour les rayonnements reerais par 
la Terre, valable tout au moins pour la region qui avait ete survolee par ce 
satellite, c'est-a-dire comprise entre les latitudes 30°N et 30°S. 

Par comparaison avec les resultats precedents (voir Fig. 1), nous avons ete ame- 
nds a decrire le processus de reemission de la fapon suivante : 

Appelons I la constante solaire (1360 W/m* lorsque la distance terre-soleil est 

egale a 1 U.A.), E. l'energie solaire incidente sur un element de surface situe 

inc 

au sommet de l'atmosphere terrestre, et la distance zenithale du soleil par 
rapport a cet element (voir Fig. 2) 


Une partie de cette 

sera reemise (E ) 
rein 


E. = I x cos 
inc 

= 0 

energie incidente 
Si l'on designe 


K si £ > 0 

0 ^O 

si £ <0 

9 

sera absorbee (E , ) 
abs 

par p 1' albedo 


( 1 . 1 ) 

tandis que 1' autre 


'abs 


= (1 - p) x E. 


inc 


( 1 . 2 ) 


E = p x E . 
rem inc 


(1.3) 


Par ailleurs, la terre emet son propre rayonnement infra-rouge E . Si l'on s'in- 

IK 

teresse au bilan radiatif terrestre, on ecrira pour un element de surface donne 


E bil 


E t - E T „ 
abs IR 


(1.4) 


L'energie infra-rouge rayonnde par la terre et l'albddo dependent de la latitude 
du lieu considdre et de la saison. On dcrira : 


'IR 

P 


f (latitude, date) 
g (latitude, date) 


(1.5) 

( 1 . 6 ) 
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On peut estiraer ces deux quantites a partir de tables. Nous avons utilise celles 
publiees par J.S. ELLIS et T.H. VONDER HAAR, 1976. 


II nous est apparu necessaire de considerer le rayonnement reemis par la terre 
comme ayant deux composantes : 


- une composante liee a un processus de rediffusion suivant la loi de LAMBERT. 

C'est l'energie rediffusee E . , 

red 

- une composante liee a la reflection speculaire tiu rayonnement solaire sur la 

surface terrestre. C'est l'energie reflechie speculairement E . 

Kb 

La difficulte est de fixer la part relative d'une composante par rapport a 1 'au- 
tre. Si on appelle a le pourcentage d'energie reflechie speculairement, tel que 


J RS 


a x E 


rem 


et done 


E = (1 - a) x E 
red rem 


(1.7) 

( 1 . 8 ) 


nous avons trouve pour la region tropicale, a partir des donnees de CACTUS, une 
representation satisfaisante de l'effet observe en ecrivant 


o = 0.1 x (1 - cos £ 0 ) (1.9) 

On verra que la valeur choisie pour o est cruciale et qua sa modification even- 
tuelle en fonction de la latitude peut changer de facon importante les resultats. 

Appelons ?«.la distance zenithale du satellite par rapport a l'element de surface 
considere 

et d la distance de cet element de surface au satellite. 

Si l'on admet la loi de LAMBERT, le flux rediffuse en direction du satellite 
par un element de surface terrestre donne dx s'exprimera sous la forme 

d<J> = p x E. x (1 - o) x cos <P x — x dT (1.10) 

red me s it d J 

Pour connaltre le vecteur 4> re( j a un instant donne, il faudra integrer ce flux 
elementaire sur tous les elements de surface de la calotte terrestre eclair^e. 

Le calcul du flux 4> refl^chi speculairement et qui atteint le satellite est un 
simple probletne geometrique. II n'y a qu'un seul point de la surface terrestre 
qui emet ce flux en direction du satellite. 
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De meme, le calcul du flux infra-rouge terrestre <J> frappant le satellite ne 

IK 

presente aucune difficulte et est tres rapide puisqu'on peut considerer ce der- 
nier corame radial. 

2. Calcul des perturbations de l'orbite sous l'effet des rayonnements infra-rouge 
et reemis 


Pour mettre ces perturbations en evidence, nous avons calcule par integration 
numerique une premiere orbite pour un satellite soumis uniquement a la pression 
solaire directe pendant une periode de 2 jours, puis une deuxieme orbite pour un 
satellite soumis en plus a la pression de radiation terrestre infra-rouge et a 
la pression de radiation reemise. Par comparaison de ces deux orbites nous en 
avons deduit devolution des elements, c'est-a-dire da/dt, de/dt, etc. 

Cette integration numerique se fait en utilisant les 6 equations de GAUSS. Nous 
rappellerons la premiere et la troisieme equations : 


da 

dt 


2 



[ Re sin v + (1+e cos v) S] 


( 2 . 1 ) 


di = r cos (oo+y) (2.2) 

dt n a 1 'J 1 -e 2 


R etant la composante radiale de la force perturbatrice 

S sa composante perpendiculaire au rayon vecteur dans le plan de l'orbite 
W sa composante perpendiculaire au plan de l'orbite. 


Pour proceder a cette integration il suffit done, a chaque pas, de calculer les 

composantes R, S et W relatives au flux rediffuse <£ re( j* au flux reflechi specu- 

lairement <f> „, et au flux infra-rouge terrestre <j> . C'est evidemment le calcul 

Kb IK 

de 9 re{ j qui ralentira considerablement le programme. 


3. Accroissement de la rapidite de calcul 

Nous nous sommes interesses au cas des satellites presentant une excentricite de 
leur orbite suffisamment faible pour que le rapport entre l'altitude de leur 
apogee et celle de leur perigee soit inferieur a 1.3. 

Si l'on veut proceder a un calcul d'orbite s'etendant sur plusieurs jours, il 
est avantageux de 1 'organiser de la facon suivante : 





Appelons 6^ la declinaison du soleil. On verifie que, sur toute la periode d'in- 
( tegration envisagee, 6 q ne varie pas plus de 5° ; sinon on decoupe cette periode 

en un nombre suffisant de sous-periodes . 

Appelons t l'heure solaire du point subsatellite a un instant donne. 

Appelons <P la latitude du point subsatellite a un instant donne. 

Les tables d' albedo que nous avons utilisees ne tiennent pas compte de 1^ lon- 
gitude des elements de surface terrestre intervenant dans le calcul. La force 
perturbatrice engendree par l'energie solaire reemise ne dependra pratiqueraent 
que des positions relatives Soleil-Terre-Satellite , c'est-a-dire du couple (t, <p) . 
II y a quasi-stationnarite de ces effets tant que la declinaison du soleil 6 q 
variera suf f isaxranent peu (on admet une variation de cette declinaison de + 2“5 
pendant la sous-periode consideree). Par consequent, il n'est pas necessaire de 
calculer cette force perturbatrice a chaque pas. On peut le faire en debut de 
programme en supposant que la declinaison du soleil reste constante pendant toute 
la sous-periode consideree, et egale a la declinaison moyenne 6^. 

a) Nous allons initialiser 6 tableaux a 2 dimensions. 

( Les trois premiers representent les valeurs des composantes R (t, <P ) , S (t, <p ) 

et W (t, <P) , pour l'altitude moyenne ho du satellite, de la force perturba- 
trice subie par celui-ci a cause des rayonnements rediffuse et speculaire. 

Les trois suivants representent les valeurs, a cette meme altitude ho, des deri- 
vees de ces composantes en fonction de l'altitude, c'est-a-dire 

— ( t <p) ££ ( t <p) — (t <P ) 

dh 1 dh V dh 9 1 

Ces six tableaux sont obtenus en faisant decrire a un satellite polaire fictif 
et pour deux altitudes differentes les couples (t, <p ) . 

t variera de 0 a 12 heures par pas de 20 minutes, ce qui represente 37 valeurs. 
(II est inutile de faire varier t de 12 a 24 heures, car, par raison de symetrie, 
on retrouve les meraes valeurs pour R et S, au signe pres pour W) . 

^variera entre -i et + i (i etant 1 ' inclinaison de l'orbite du satellite etudie), 
avec un pas de 5°. Ceci nous conduit a un maximum de 37 valeurs pour un satellite 
polaire. 


I 
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b) II est maintenant possible de calculer a chaque instant et dans le systeme de 

reperes arbitraire des tableaux la force subie par le satellite par simple 

interpolation lineaire, apres avoir corrige les coordonnees du satellite (t , <P ) 

en fonction de la declinaison exacte du soleil 6 . Par les rotations necessai- 

e 

res, nous obtenons .1 'expression du vecteur acceleration du au rayonnement reemis 
pour le satellite considere. 

En ce qui concerne <f> , nous procedons au calcul direct a chaque pas. C'est un 

XR 

calcul tres rapide. 

Grace a cette organisation, nous avons obtenu, avec un pas d ' integration de 1 mn 
et pour un satellite polaire, une rapidite de calcul multipliee en moyenne par 10 
avec une difference entre les resultats obtenus pas a pas et ceux obtenus par 
cette methode de l'ordre de 0,2Z. Le gain en rapidite est encore augmente si 1* 
inclinaison de l'orbite est plus faible. 

4. Methode analytique 

Dans le cas d'un satellite ayant une orbite evoluant lentement (satellite geode- 
sique a haute altitude par exemple) , il est possible, si les perturbations dues 
aux rayonnements reemis et infra-rouge ont ete calcul£es sur une periode suffi- 
samraent longue, de faire une analyse de FOURIER de ces perturbations. A partir 
de cette analyse, il sera possible d'estimer pour toute date demandee et avec 
une faible imprecision, la valeur de la perturbation par la somme de quelques 
termes sinusoldaux multiplies par des constantes. Cette derniere procedure, qui 
necessite d'abord 1 'application de la methode precedente, est de loin la plus 
rapide, d'ou son interet. 

5. Application a LAGEOS 

1) A partir du calcul des perturbations de l'orbite de LAGEOS sous l'influence 
des rayonnements reemis et infra-rouge pour 60 dates espacees d'un mois, 
grace au procede decrit au paragraphe 3., nous avons fait 1' analyse de FOU- 
RIER des resultats. 

En ce qui concerne le demi-grand axe, nous avons trouve que son evolution pouvait 
etre decrite par 
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( 


— = [ 0.063 sin (2A - ft + 10°) 

dt 1 e 

+ 0.049 sin (2A - 2 ft + 6°) 

© 

+ 0.029 sin (A - ft - 97°) 

+ 0.053 sin (ft + 10°) ] ram/jour 

A etant la longitude du soleil 

o 

et ft 1' ascension droite du noeud ascendant du satellite. 


2) Les valeurs observees de 1 'acceleration de LAGEOS le long de sa trajectoire 
ont ete determinees par plusieurs auteurs (SCHUTZ B.E., TAPLEY B.D., 1980, 
SMITH D., DUNN P. , 1980). 


La comparaison avec les valeurs de la derivee du demi-grand axe — est possible: 
la premiere equation de GAUSS (2.1) permet d'ecrire dans le cas de LAGEOS, puisque 
e = 4. 10" J , 



da 

dt 


Si l'on examine les resultats du calcul, on remarque que l'on ne retrouve ni les 
amplitudes, ni la forme de la courbe des valeurs observees (voir Fig. 3). 

A la suite de cette constatation, nous avons refait les calculs en augmentant 
1* importance du rayonnement reflechi speculairement au detriment du rayonnement 
rediffuse (Eqs 1.7 et 1.8). L'amplitude des perturbations calculees pour le demi- 
grand axe de LAGEOS a augmente presque proportionellement a la valeur atteinte 
par le coefficient a aux poles.,. II est possible de retrouver l'ordre de grandeur 
des perturbations observees si l'on admet un coefficient o proche de l'unite, ce 
qui montre le role important que pourrait avoir le rayonnement reflechi dans la 
perturbation de l'orbite de LAGEOS, ainsi que l'ont signale ANSELMO L. et al., 
1983. Toutefois, la forme de la courbe n'a pas dvolue sensiblement vers celle 
des valeurs ■ observees . 


6. Conclusions 

Nous avons calcule les accelerations perturbatrices dues aux rayonnements reemis 
et infra-rouge que subirait le satellite LAGEOS sous 1' influence d'un modele de 
rayonnement construit a partir des accelerations correspondantes mesurees a bord 
du satellite CASTOR dans la region tropicale. Les resultats des calculs sont fort 
differents des valeurs observees. 
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Si nous supposons que les accelerations observees pour LAGEOS sont dues aux 
rayonnements reemis et infra-rouge, il faudrait admettre qu'un modele de rayon- 
nement valable pour CASTOR en zone equatoriale ne le serait pas pour LAGEOS qui 
survoLe aussi Les zones polaires. Nous avons pu constater, grace a des simula- 
tions, 1 ' importance que joue le rayonnement reflechi speculairement pour des 
satellites ayant une inclinaison elevee d'orbite. On reraarque aussi qu'un modele 
d'albedo symetrique (par exemple p = a + b x sin 2 <p ) est trop simpliste dans ce 
cas, car les forces perturbatrices calculees sur une orbite tendent a se corapen- 
ser . 

Peut-etre pourrait-on trouver un meilleur accord entre les resultats du calcul 
et les valeurs observees en supposant que le coefficient de specularite 0 varie 
dans le temps, par exemple en fonction de la saison locale ou de la raeteorologie . 
II est certain que si l'on admet que ce coefficient n'ait pas la meme valeur 
a un instant donne dans les zones polaires arctique et antarctique, on peut abou- 
tir a des accelerations le long de l'orbite encore plus importantes. Cette hypo- 
these semble difficile a tester avec les donnees existant actuellement . Mais nous 
touchons la un point crucial. 

Peut-etre faudrait-il plutot chercher aux perturbations de l'orbite de LAGEOS 
une autre origine principale. Nous savons deja que la decroissance seculaire du 
demi-grand axe de son orbite est tres vraisemblablement liee a la charge electri- 
que du satellite, comme cela a ete propose par plusieurs auteurs (AFONSO G.B. et 
al. 1980, RUBINCAM D.P. 1980). On peut penser que les parametres electriques 
subissent des variations periodiques. Mais la connaissance insuffisante de ces 
parametres ne permet pas encore actuellement de chiffrer avec certitude leurs 
effets. 

Par ce travail nous avons reussi a augmenter dans une proportion interessante la 
rapidite de calcul des effets du rayonnement reemis. Nous avons montre le role 
important que joue le rayonnement reflechi speculairement pour les satellites 
ayant une inclinaison d'orbite elevee. S'il est encore impossible d'etablir 
une solution definitive, nous avons propose un modele ''classique" de rayonnement 
reemis. II sera tres facile d'adapter le programme correspondant des qu'une con- 
firmation des valeurs du coefficient de specularite O sera obtenue. 

Des etudes et des analyses de donnees etalees sur une plus grande periode devront 
etre poursuivies, afin de determiner dans quelles proportions chacune des deux 
hypotheses precedentes interviennent pour expliquer les observations. La suite 
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normale de ce travail sera de tenter une representation coherente du phenomene 
observe par la juxtaposition d'un modele de rayonneraent reemis et d'un modele 
de freinage electrique. 
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